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Seguir-se-á a descrição teórica da associação em paralelo de dois rectificadores com 
retorno pelo neutro trifásicos, alimentados por uma dupla estrela trifaásica. A generali- 
zação dos resultados obtidos permite a análise e o dimensionamento da montagem 
apresentada na fig. 1. 


2.1 Rectificador de Dupla Estrela Trifásica 


A montagem analisada é o rectificador de dupla estrela trifásica, constituída pela as- 
sociação em paralelo de dois rectificadores com retorno pelo neutro trifásicos de cáto- 
dos comuns alimentados por um transformador de dois secundários trifásicos cujos 
enrolamentos são ligados em oposição de fase, figura 2. Nesta figura, o enrolamento 
do primário está ligado em estrela, no entanto a ligação triângulo pode ser também 
utilizada, como se mostra na figura 1. 


Se os rectificadores são comandados, então os ângulos de atraso à condução dos tiris- 
tores (£) devem ser idênticos de modo a estabelecerem tensões de saída com o 
mesmo valor médio. No entanto as fontes que alimentam os dois rectificadores não 
estão em fase, e assim é necessário colocar entre os dois rectificadores uma indu- 
tância de ponto médio - bobina interfases - aos terminais da qual é aplicada a dife- 
rença instantânea das tensões de saída dos dois rectificadores. 


A análise desta montagem que a seguir se apresenta, é baseada na referência [3]. 


Em funcionamento normal, figuras 3a e 3b, a tensão de saída é igual a: 


o Up MU 
a = = (1) 


figura 2 
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A tensão u, aos terminais da bobina interfases é dada por: 


Deste modo, a tensão de saída da montagem, uy, é dada por: 


u, + us 
RU 
Es1] 
U1 
=“ % 
Fa PA, »* se Pi * 
/ * “ * 1 f 
z P a e 4 = Me e “4 
% * a £ % 
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* a % E. "4 * 
poi E nao? O | ndo RARE À PR > 
1 Vosnh 
Esi2 o 
Us 
4% Ci % E 
4 Pd A ie sr CY us Tm Sho 
" £ 4 ,” er TR 
*w4f " dá ko Vá Pd “4 Fá 
Pao , rs ns. a. 
= a a ! rd % 
n.º af PAD du DO qi Cm a” OR Já a 


figura 3a 
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Desprezando as reactâncias de comutação, o comportamento da montagem pode ser 
descrito através das funções de existência ou de comutação feptool) e fpjl(ot) com 
je [1,2,3] que caracterizam respectivamente o estado dos tiristores dos conversores 1 e 
2. Elas podem tomar os valores O ou 1; o valor 1 corresponde ao estado condutor do 
tiristor respectivo, enquanto que o valor O significa bloqueio. 


Considerando que as tensões das fases do primário são dadas por: 


| po ATE 
tu =Vam cos 0x =2j 5] (4) 
| 4 
podemos escrever: 
= vd : = GR 
fmtor)=16512j-3] 7 +e+2kT Sm <S(2j-1+E+2KT (5) 
| 3 
E: : ro roÃá À 
fotnj=1652;-3) q TEHÍZK +lir<ws<(2;-1) q FE(ZK +lym (6) 


com K=-0. 1,2, 3: 


As tensões de saida dos dois conversores são dadas por: 


lj = 3 Fatores to (7) 


ly = Y ir tor es, or) (8) 
= 


Tendo em conta que 
contr) = ec 0%) 


a tensão de saida da associação em paralelo é dada por: 
ly => RO p= fem OM esto (9) 
aa = 


A função fp Áwt) é idêntica à definida para a ponte de rectificação trifásica com fonte 
em estrela. Assim, este último conversor, à parte o factor 1/2, apresenta uma tensão 
de saída idêntica à da montagem em análise. Em particular, a tensão rectificada apre- 
senta um indice de pulsação p-6. 


A corrente de carga (!,) divide-se em duas, cada uma destas associada a um rectifi- 
cador. Uma corrente alternada 1, que resulta da aplicação da tensão (u,-u,) à indutân- 
cia L, soma-se e subtrai-se à corrente de cada um dos rectificadores: 


MP era 
== + (10) 
a go a 
= a (11) 


Como as correntes 1, e |, não podem ser negativas, o funcionamento que se descre- 
veu so é possivel se o valor da corrente de carga é superior ao dobro do valor máxi- 
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mo da corrente i,. O valor da corrente ly abaixo do qual não se pode descer em funci- 
onamento normal, é chamado o valor da corrente crítica. 


E de assinalar que as componentes continuas iguais a //2 que circulam nas duas 
metades da indutância não provocam nesta uma magnetização continua, elas circu- 
lam em sentidos opostos de modo que a sua contribuição para a força magnetomotriz 
é nula. 

O valor médio das correntes secundárias associadas a cada uma das fases do trans- 
formador é nulo: 


E mo sam = Q (12) 
Ho 


Desprezando a corrente de circulação que na prática apresenta um pequeno valor face 
ao máximo da corrente que a montagem pode debitar, as correntes do primário são 
idênticas às obtidas com uma montagem trifásica em ponte. 


2.1.1 Análise da tensão de saida 


A análise das figuras 3a e 3b mostra-nos que apesar das tensões u, e u, terem o 
mesmo valor médio, instantaneamente elas são distintas, u,-U5*0. Este facto estabe- 
lece uma corrente de circulação i,. Esta é determinada, fundamentalmente, pelas ter- 
ceiras harmónicas presentes à saída de cada um dois rectificadores com retorno pelo 
neutro trifásicos. 


Desprezando a queda de tensão média na bobina, o valor médio da tensão de saida 
da associação em paralelo é igual ao valor médio da tensão de saida de cada um dos 
rectificadores que a constituem: 


— cos E) (13) 


As componentes harmónicas de ordem 3 (e multiplas) geradas pelos rectificadores 
com retorno pelo neutro trifásicos não aparecem à saída da sua associação em parale- 
lo. Esta caracteriza-se por apresentar um índice de pulsação p=6, isto é, a frequência 
angular da harmónica de ordem mais baixa é 60 (300 Hz). 


Do ponto de vista da tensão de saida, a associação em paralelo de dois rectificadores 
com retorno pelo neutro trifásicos, representada na figura 2, é equivalente a um rectifi- 
cador com retorno pelo neutro hexafásico, isto é, de seis fases (n=p=6), cujas forças 
electromotrizes têm uma amplitude Vy equiv. 1º 


Com esta equivalência a tensão de saida da montagem é dada: 


ty ae Hy, + pa Hj (1 4) 
b=| 


onde a componente contínua, isto é, o valor médio é dado por: 


n[7) 
sin) 
6 


Hym = V We esqurv r COSsE (1 5) 
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a componente harmônica de ordem k6 é dada pela expressao: 


His = Cr cost kom =D) (16) 
com 
sin] Es | 
2V E: | | o o 4 " 
4 anda A fi S 6, Vlk6) sin“e+cos e (17) 
Ik6/ —1 | 
6 


A amplitude equivalente, Vyequiv 1 é Obtida tendo em conta a igualdade dos valores 
médios das diferentes tensões, (13), e as suas respectivas expressões, (13) e (15). 
Estes resultados estabelecem: 


WE 
sn sn 
6 ? 
Vau e tose=V.. ---Stose 18 
HW equiv Tt MH T ( ) 
6 3 


figura 4 


o que conduz a: 


| 


Wequivi — T | (19) 
2.1.2 Análise da corrente no primário do transformador 


A análise das figura 3a e 3b mostra-nos que as correntes nas fases do primário do 
transformador têm uma forma semelhante à estabelecida pelo rectificador trifásico em 
ponte. Desprezando os efeitos da corrente de circulação, a associação em paralelo de 
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dois rectificadores com retorno pelo neutro trifásicos apresenta no primário do trans- 
formador uma corrente constituída por duas ondas rectangulares de largura 27/83, 
desfasadas entre si de mr, uma de amplitude (+//2) e a outra de amplitude (-1/2). 
Estas ondas estão desfasadas de e relativamente às tensões primárias corresponden- 
tes e cada uma está associada a um rectificador com retorno pelo neutro. 


As harmónicas presentes na corrente de fase no primário são de ordem impar (2k-1): 


à ox )=— sim mt-e-2j-—l 


My; 3 4 - | [ (2% =D | asl(Di E | E |] (20) 
nem 3 4” L 33) 


com ke [1,2,3, . Je je [1,2,3). 


À parte, a fundamental de pulsação igual à da rede (50 Hz), a harmónica de ordem 
mais baixa tem uma frequência angular igual a 5 (250 Hz). 


3. Montagem com transformador em triângulo 


Como já se fez referência e a figura 1 ilustra, o transformador com secundário hexa- 
fásico que alimenta os dois rectificadores com retorno pelo neutro trifásicos pode ter o 
primário ligado em triângulo. O comportamento da montagem resultante é semelhante 
ao estabelecido pela associação de dois rectificadores com retorno pelo neutro trifasi- 
cos com transformador cujo primário está ligado em estrela, figura 2. 


Os comportamentos das duas montagens, respectivamente com primário ligado em 
estrela e com primário ligado em triângulo, alimentadas pela mesma rede, apresentam 
diferenças: 


as tensões de saída das duas montagens têm componentes flutuantes desfasadas 
entre si de 180”, como resultado da desfasagem entre os sistemas hexafásicos à 
entrada dos rectificadores propriamente ditos; 


- as correntes pedidas a rede, tendo o mesmo conteúdo harmónico, apresentam no 
entanto formas de onda distintas: uma é corrente de fase (estrela) e outra resulta 
da diferença de duas correntes de fase (triângulo). 


4. Montagem dodecafásica 


Com a metodologia apresentada no parágrafo 2 mostra-se que o esquema represen- 
tado na figura 1 é uma montagem dodecafásica. Assim, o índice de pulsação da ten- 
são de saída é p=12 e a corrente pedida à rede, teoricamente, não apresenta harmó- 
nicas de ordem inferior a 11. Na figura 5 ilustra-se o comportamento da montagem 
dodecafásica. 

Na ligação em paralelo de duas associações de rectificadores com retorno pelo neutro 
trifásicos, uma com o primário do transformador em estrela (Y) e a outra com o primá- 
rio do transformador em triângulo (A), as tensões de saída de cada uma devem apre- 
sentar o mesmo valor médio: uy,=Usm: Este facto é conseguido se as amplitudes das 
tensões dos dois sistemas hexafásicos forem iguais e se o ângulo de comando for 
idêntico para todos os tiristores. 


Entre os dois sistemas hexafásicos existe uma desfasagem de 30º, deste modo ape- 
sar das tensões uy e u, terem o mesmo valor médio, instantaneamente elas são 
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distintas, uyzu,. Este facto estabelece uma corrente de circulação i,; de frequência 
angular 6 (300 Hz), cujo valor é definido pela bobina interfases, L,. 


qm 


qm 


Am 


qm 


figura 5 


Desprezando a queda de tensão média na bobina L,, o valor médio da tensão de sa- 
ida da montagem dodecafásica é igual ao valor médio da tensão de saída da cada 
uma das duas associações em paralelo de rectificadores com retorno pelo neutro tri- 
fásicos que a constituem: 


UH, = Hy = EE am (21) 


As componentes harmônicas de ordem 6 (e múltiplas) geradas pelas associações em 
paralelo de rectificadores com retorno pelo neutro trifásicos não aparecem a saida da 
montagem dodecafásica. Do ponto de vista da tensão de saida, a associação em para- 
lelo de duas associações em paralelo de dois rectificadores com retorno pelo neutro 


ES == 126 


E Em Eis EINE E Sa SR 


revista de Engenharia - 92 


João Jose Esteves Santana 


trifásicos, representada na figura 1, é equivalente a um rectificador com retorno pelo 
neutro dodecafásico, isto é, de doze fases (n=p=12), cujas forças electromotrizes 
têm uma amplitude Vy equiv.2- Com esta equivalência a tensão de saída da monta- 
gem é dada: 


U=u, + * Ui (22) 
h=] 


onde a componente continua, isto é, o valor médio é dado: 


sin 
12 


UH, — Vy equir? o cose (23) 
com 12 

ME 

sin[ 3 ] 
Vu equiv? e Vu fr (24) 

4sin[ ] 

12 

Us = Cu cos(k120r — Ou) (25) 


com 


sin| 7 
2V u equivo | 2 


(k12) -1 T 
12 


Em cada uma das fases da rede, a corrente é o resultado da soma de duas correntes: 
uma associada ao primário do transformador em estrela e a outra é a corrente de linha 
do primário ligado em triângulo. Deste modo há harmónicas que se anulam; à parte a 
fundamental de pulsação angular w (50 Hz), a harmónica de mais baixa ordem tem 
uma pulsação igual a 110 (550Hz). As correntes nas três fases da rede estão desfasa- 
das entre si de (27/3); a corrente de rede na fase 1 é dada por: 


Can = (26) 


2 Pe IR di NR A PN am BECO o E 
a 1] sin) , [costa — E) + Sin : Jeos[t lo -e)]+ sin : Jeosft% or + e (27) 


/ 
nr é | 5 | 3 NR: 


Lirgundr 


Libano Monteiro (FRAEP) e José Carlos Lameira (PRINEC) reconhece-se o 
profissionalismo demonstrado na realização da fonte e no projecto da ligação à rede 
da EDP, respectivamente. 
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5. Conclusões 


Neste artigo foi analisado o circuito electrónico de potência da fonte de alimentação 
das bobinas que criam o campo magnético toroidal do TOKAMAK do Centro de Fusão 
Nuclear do IST. 


Atendendo às caracteristicas exigidas à saida (baixas tensões e elevadas correntes), 
foi seleccionada a associação em paralelo de rectificadores com retorno pelo neutro. 


Desprezando as reactâncias de comutação, a definição de funções de existência, ca- 
racterizando o estado dos semicondutores, permitiu a elaboração de um programa de 
simulação da montagem. Este e o desenvolvimento analítico conduzem à caracteriza- 
ção do funcionamento e dimensionamento das diferentes componentes do sistema. 


Actualmente a fonte encontra-se em operação, e assim futuros artigos deverão ser 
apresentados, nomeadamente os relacionados com a interligação do circuito electró- 
nico de potência à rede da EDP e ao TOKAMAK, os sistemas de protecção e controlo 
do conjunto e a caracterização experimental. 
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Resumo 

“análise por difusão de Thomson é uma 
écnica muito usada na medição da den- 
sidade e da temperatura electrónica de 
um plasma laboratorial. Neste trabalho 
apresentamos os princípios básicos & 0 
projecto do diagnóstico de difusão de 
Thomson do 'tokamak' ISTTOK. 


Abstract 

Thomson scattering of light is a com- 
monly used technique for measuring 
electron temperature and density of a 
laboratory plasma. In this paper we 
present the basic principles and the 
eneral design of the Thomson scattering 


diagnostic for the tokamak ISTTOK 
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THE THOMSON SCATTERING 
DIAGNOSTIC FOR THE 
TOKAMAK ISTTOK 


| Introduction 


Laser beam diagnostics are very often used as a means for active probing of plasmas, 
without distortion of the plasma parameters. Depending on the specific laser-plasma 
interaction involved (absorption, refraction or scattering) information on a wide variety 
of plasma parameters can be obtained with high spatial and temporal resolutions. 


In this paper we present the basic principles and the general design of the Thomson 
scattering diagnostic used in the tokamak ISTTOK [1] for electron temperature and 
plasma density measurements. 


ISTTOK is a small tokamak that was assembled from the basic components of the 
tokamak TORTUR, which was discommissioned in 1988 by the EURATOM/FOM 
Association, in Niguwegein, The Netherlands. The diagnostics, the control units and the 
data acquisition system for ISTTOK have been designed and installed by portuguese 
physicists, engineers and technicians. 


2 Basic Principles 
When an electromagnetic wave interacts with a moving charged particle, a scattered 
wave is emitted showing a Doppler shift of the wavelength due to the motion of the 


charge [2]. For a large collection of scattering centres, we have to consider both space 
and velocity distributions of these charges. 


It has been demonstrated [3] that if the wavelength of the incident radiation (A;) is much 
smaller than the Debye length (Ap), so that «=A;/An<<1, the collection of scattering 
centres can be treated as having a random spatial distribution. Under these conditions, 
scattering is said to be incoherent and the scattered power has a spectral distribution 
similar to the velocity distribution of the scattering centres. Considering relativistic cor- 
rections, it can be shown [3,4,5] that the scattered power is written as 


Po SA Print | «AA “AA CÃA IA 
P(A, Jd dO = 1-3 EAR, 76% [5 exp] E Cedo (1) 
2yn sin| — | Ei 4vaÃ, sin” - | em Avi A sin” | Mg 
| É. | | si, , 2 | ei | ] 


where v, is the thermal velocity (2kT/m,)!2, AA is the wavelength shift, P, is the inci- 
dent power, rp is the electron's classic radius (15=2.8196x10!9 cm), L is the length of 
the scattering region, n, is the plasma density and O is the scattering angle. The expo- 
nential term is a consequence of assuming a maxwellian velocity distribution. In this 
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For the low electron temperatures expected in ISTTOK no relativistic effects need to 
be taken into account. The wavelength shift is sufficiently small to neglect the terms of 
order higher than zero in the polynomial factor. Introducing a normalised wavelength 
shift 


g=1— (2) 


we have 


P dx dQ = id exp E 2, um dv dO (3) 


which has a convenient form to show, except for a scale factor, the independence of P, 
irom the incident wavelength. 


The measurement of the electron temperature is based on the Doppler broadening of 
the incident wavelength. From the previous expression we have for the width at 1/e of 


the peak 
PA a sin| : g 
-2—SE 7911x10º VT; sin[5 ] (4) 


l 
with the electron temperature (T,) in ev. 


The plasma density can be calculated from the amplitude of the spectrum only after 
absolute calibration. Then the scattered power is retrieved from the detection signal to 
give for x=0: 


dê n L 
Jima q = 
vy 2NT sn(S) JO Ela 


Sn 


where np and L are expressed in cgs units and Ta in ev. 


3 Expected range of plasma parameters 


A portuguese team made about sixty shots in Nieuwegein, operating TORTUR at 
conditions similar to those which are expected for the ISTTOK operation: a toroidal 
magnetic field equal to 0.45 T and a plasma current of about 5 kA. 

Figure 1 presents a typical output obtained with the TORTUR Thomson scattering di- 
agnostic [6]. 


Figure 2 shows the histograms of the measurements (T, and no) made with this diag- 
nostic. From these data we expect that, also in ISTTOK, the electron temperature shall 
be in the range 


1 STS LIDEV 
and the plasma density in the range 


IxiQt cms PAO gm 
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figure 1 Atypical Thomson scattering output 
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figure 2 Histograms of the electron temperature and the plasma density measurements made with the Thomson scattering diagnostic 
of the tokamok TORTUR. 


For these ranges of plasma parameters, the condition for incoherent scattering is ful- 
filled, since a=A/Ap ranges between 10% and 102. For the low temperatures involved 


the wavelength spectrum of the scattered light can be represented by a gaussian with- 
out relativistic correction. 
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4 General set-up for the diagnostic 


In this section our attention is focused on the general layout of this diagnostic, the 
geometric constraints imposed by the experimental device, the spectrometer and de- 
tector design, the laser characteristics, and, finally, the data acquisition interface. 


4.1 General layout 


The experimental set-up of the Thomson scattering diagnostic for a measurement at 
one position in the plasma is presented in Figure 3. 


INPUT OPTICS 


PLASMA 


NA AT 
SE RI rs 


DATA AQUSITION 


LASER DUMP 
SPECTROMETER AND DETECTORS 


figure 3 The experimental set-up of the Thomson scattering diagnostic 


The incident electromagnetic power is provided by a laser beam, which is guided by a 
series of suitable mirrors and crosses the plasma column. Generally, ruby and neo- 
dynium lasers are used in this type of diagnostic. For our purpose a Nd:glass laser 
(1054 nm) was chosen to avoid its output frequency overlapping the spectral range of 
impurity line radiation from the plasma (the radiative transitions are mainly in the 
neighbourhood of the ruby laser frequency (1=694.3 nm)). 


The light emitted by the scattering volume is collected within the solid angle by a lens 
system and guided into the spectrometer, which is connected to the data acquisition 
system. The length of the scattering region (L) is limited by the aperture of the optical 
elements and determines the spatial resolution. Generally, the lens system is chosen 
such that the scattering volume is imaged on the detector dimensions. 


For ISTTOK a scattering angle of 90º and L=1 cm were chosen. 


We now evaluate the maximum solid angle A£) allowed by the tokamak ports. There 
are three ports available in ISTTOK for Thomson scattering measurements (figure 4), 
with a diameter of 38 mm. Window 1 can be used for profile measurements. It has the 
largest exit port diameter of the three and the laser beam crosses the plasma at a mi- 
nor radius of 22.5 mm. This allows a high solid angle at a small distance from the 
plasma centre. Windows 2 and 3 have a smaller exit port diameter. In the former the 
laser beam crosses the plasma at 1=60 mm, in the latter at 5 mm. The maximum 
solid angles for each viewing port are given in the Table 1. These values will be used 
for calculations and system dimensioning. 
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figure 4 — Schematic view of the ISTTOK diagnostics 


Table 1 
| Por o r(mm) “r(mm) | AQ (mst) 
| 225 | 236 
| 
| 


60 
3 


In the first phase of operation of this diagnostic we will perform single point measure- 
ments. By doing measurements of successive and reproducible discharges, obtained 
with the same external parameters, we expect to obtain profiles of the plasma density 
and electron temperature, alternately using the three entrance poris at 1=5, 22.5 and 
60 mm. 


In the second phase we will do measurements simultaneously at two points in the same 
discharge, using a beam splitter to distribute the laser over the two entrance poris at 5 
and 60 mm. In this phase two collecting systems and two spectrometers are needed. 


Finally, in the third phase we will do measurements at four points (at two different ra- 
dial and vertical positions) in the same discharge, again using a beam splitter and two 
entrance ports. However in this case four collecting systems and two spectrometers, 
each one with two inputs, are needed. The two collecting systems are connected to the 
same spectrometer through fiber bundles with different lengths in order to create an 
adequate delay between the two signals. 


Since the minimum time between two consecutive laser shots is much greater than the 
duration of the ISTTOK discharges, we can only do single time measurements. The 
temporal variations of the plasma density and electron temperature will be obtained 
from successive and reproducible discharges by varying the instant of time at which 
the laser is fired. 
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4.2 Optics 


There are two optical systems in this diagnostic. The first one is to transport the laser 
beam to the entrance port, to image the beam onto the scattering volume and to dump 
the laser beam after crossing the plasma column; the second one is for collecting and 
resolving the scattered power. 


The first system is mainly composed by mirrors (to guide the beam to the entrance 
port), focusing lens and the entering port window. In the bottom side of the vacuum 
chamber there is a laser dump. Special diaphragms are used in both sides to avoid the 
stray light. These optical components must be adequate to sustain the very high in- 
tensity of the laser beam, which is necessary regarding the low densities expected in a 
common ISTTOK plasma. 


The second system includes the output port window, an objective lens, a fiber bundle, 
a spectrometer and the final lens in front of the detectors. A similar study [7] showed 
that it is reasonable to expect an overall transmission of 10 to 20%. 


4.3 Spectrometer 


As was previously referred (equation (5)), the plasma density is, roughly speaking, 
calculated by measuring the amplitude of the scattered light, and the electron tempera- 
ture from the spectral width. We have, therefore, to measure the scattered power 
transmitted to a number of spectral channels in order to reconstruct the spectrum. A 
trade-off does exist between spectral accuracy (high number of channels), feasibility 
(since many channels imply a low signal level in each one) and price. 


Assuming that the plasma is in thermal equilibrium, the spectrum of the scattered light 
must be a gaussian (equation (3)) that can be obtained from the amplitudes given by, 
at least, two channels [8]. A third channel is also used to provide confirmation of the 
assumption of the plasma thermal equilibrium. 


Since we are not absolutely sure to have always an equilibrium distribution, it is 
thought useful to measure more accurately the scattered spectrum, specially in some 
critical cases, like heating and current drive with electron cyclotron waves. In this case 
the distribution function is likely to be distorted. For that, we need at least 20 spectral 
channels. To cover the range 50 < T, < 300 eV and for the Nd:glass laser wavelength, 
the required spectrometer should be able to measure the range 1054 + 150 nm, with a 
resolution better than 2 nm for each channel. 


Spectrometers, generally, consist of a diffraction grating or a sequence of interference 
filters, and detectors. In the first case the light is dispersed according to its wavelength 
and transmitted to an appropriate number of detectors. In the second case light within 
the passband of the first filter will be transmitted to the first detector; light outside this 
passband will be reflected to the second filter and so on. 


During the ISTTOK inductive operation, we will use a spectrometer with three interfer- 
ence filters (figure 5), covering the bands 1054-1042 nm, 1042-1024 nm and 1024-964 
nm [9]. In order to minimise the effect of stray light, the laser wavelength will be 
filtered. 


For the ISTTOK non-inductive operation, with electron cyclotron current drive, we will 
increase the number of spectral channels up to twenty, for the above mentioned 
reasons. In this case we will probably use a spectrometer with a diffraction grating [6]. 


The detectors more often used in Thomson scattering diagnostics are photomultipliers 
[6] and avalanche photodiodes [7,10]. 
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At ISTTOK we use avalanche photodiodes because they have high performance in the 
infrared region. The quantum efficiency is high (up to 80%) and the active area is large 
enough to allow good optical matching of étendue with the scattering length. 


In choosing the right type of avalanche photodiode, we have to make a trade-off be- 
tween performance and an enhanced noise level, which is due to avalanche 
multiplication. 


We have chosen RCA's avalanche photodiode C30950, or a modified version of It, as 
is used in ASDEX [11,12] as well as in TORE SUPRA [13]. lts characteristics, at 
900 nm, are 3x10º V/W , a quantum efficiency (n) of (dropping to about 30% at 1.06 


FILTER 2 FILTER 3 


DETECTOR 3 


FILTER 1 


figure 5 — Schematic drawing of a spectrometer with three interference filters 


um, about the wavelength of Nd lasers) and a noise equivalent power of around 10719 
W/Hz!?2. It can be shown [8,9] that the noise arising from the avalanche multiplication 
has the effect of lowering the quantum efficiency to give an effective efficiency Nag =N 
![F, where F is the excess noise factor due to the avalanche multiplication. F varies 
from 5 to 10 for a gain ranging between 25 and 100 [7,10]. 


4.4 Signal-to-noise ratio 


Integrating equation (3), the power scattered from a plasma region with length L within 
a solid angle A£2 and transmitted into a detector can be written as : 


L 
JT. sm : ] 


where we have included the incident laser power (P) the overall transmission (7) of the 
optical system, the solid angle (442), and the spectral with (Ax,,). The factor 1/N ac- 
counis for the desired number of channels (N) whereas the factor 1/2 results from con- 
sidering a channel detecting on the average only half the maximum spectral power. 
Using the expected values of density and temperature, a laser power of 100 MW, and 
the geometric factors mentioned, we obtain typical values of L=1 nW for N=20. 


| pa 


| 
: (6) 


E = 1.13x10" Ba, T AQ Ax,. 
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The light emitted by the plasma itself (Bremssirahlung, line radiation due to electronic 
transitions and recombination radiation) will give an undesired background contribution 
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to the total detector output. We can only give estimation of the Bremsstrahlung since 
the discontinuous spectrum of bound-bound or free-bound transitions strongly depends 
on the impurities in the ISTTOK plasma, for which we have up to now no estimate. 


The Bremsstrahlung level for a hydrogen plasma is given by [3]. 

TA |.24x 10" ASS 
5 8 eXp|-—— Iv, | 
A; JT. A;T. “47 
We can make the rough assumption that the total plasma light has about 100 times the 
Bremsstrahlung level (without taking into account the line radiation from electronic 


transitions). If we choose an incident wavelength in the infrared region, we can expect 
the levels of transition radiation to be negligible. 


B, dS2 dA. —2 09x 10 


To estimate the ratio of scattered to plasma background light (signal to noise ratio 
(SNR)), we will use an expression derived in [7] and [10], which includes the scattered 
power L measured by each detector, the power level of plasma light (Po), the noise 
equivalent power (NEP) and the effective quantum efficiency (Neg): 


(8) 


In that expression T is the time during which the detector is measuring scattered power. 
This means that SNR depends on the bandwidth of the detection system. Obviously, in 
order to minimise the plasma light collected during a measurement, it is convenient to 
let T approach the laser pulse length. 
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figure 6 — Variation of the signal-to-noise ratio with the plasma density, for different laser energies. Other parameters: T=2x10? 5, 
N=20, n=0.3; |) 1=0.2, 2=100 msr, 65; 11) 7=0.2, 2=40 msr, 5; 1) 750.2, 92=40 msr, 10; IV) 1=0.1, 02=40 msr, 
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With the numerical assumptions mentioned throughout this paper, and for a Nd:glass 
laser, the SNR has a dependence on the plasma density which is shown in figure 6. 
Curves have been calculated for different laser energies (5, 10 and 15 J), and with 7 
equal to the pulse length. 


The most optimistic, still realistic, situation is presented in figure 61. It corresponds to a 
transmission factor of 20%, a solid angle of 100 msr, and a relatively low avalanche 
gain excess noise F of order 5. With these parameters, we can see that a SNR of over 
20 dB (which is an adequate value for the detection system) can be obtained with a 
10 J laser for a plasma density of about 10!“ cm. Figures 6II, 6IIl and 6I1V correspond 
to less optimistic parameters. If the transmission is as low as 10%, a solid angle of 40 
msr is chosen and if the signal is so weak that a high avalanche gain is required, imply- 
ing a high F of order 10, we will encounter the most pessimistic conditions considered. 
For them we can see (figure 6IV) that a SNR of 20 dB or higher can only be achieved 
with a 15 J laser with a detection limit of 5x10!2 cm? 


4.5 Laser characteristics 


We can also calculate the required laser energy (W) to achieve a required signal to 
noise ratio, by rewriting equation (8) 


hv | == 
SMT NEP*n, 
2:My 5 4 2hv 
hv 
AN 


As for the SNR, the required laser power dependence on some parameters described 

earlier, is shown in figures 7 e 8 represents the SNR dependence on the laser energy 
tor different plasma densities. We can conclude that detection with at least SNR=14 dB 
at n=10"2 em'3 could only be done with a 20 J laser. 
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figure 7 The required laser energy as function of the plasma density, for SNR equal to 14dB and 20dB. Other parameters: 


t=2x 10º 5, N=20, n=0.3; 1) 150.2, 2=100 msr, 5; 11) 7=0.2, 2=100 msr, F=10; HH) 7=0.2, 2=40 msr, 55; 
IV) 750.1, 42=40 msr, FO. 
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In agreement with previous considerations, we decided to buy the Quantel PG 38 
Nd:Phosphate glass laser, which has the following main characteristics: 25 J of maxti- 
mum output energy, pulse duration equal to 25 ns, 0.4 mrad of divergence at 90% en- 
ergy, beam diameter equal to 25 mm, spectral width of about 2 nm, linear polarisation 
and repetition rate of 4 shots per minute. 
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figure 8 — Variation of the signal-to-noise ratio with the laser energy for different plasma densities. Other parameters: 1=2x 10º s, 
N=20, n=0.3; 1) 1=0.1, 42=40 msr, E10; H) 7=0.2, 42=100 msr, 5. 


renas 138 


E to Es DO E Ca 


revista de Engenharia - 92 


The Thompson Scattering Diagnostic for the Tokamak ISTTOK 


o) (6) 


a 00 E o EH OC 


The schematic drawing of the PG 38 QUANTEL laser. Listing of components: 1) HeNe laser; 2) Flat total refledion mirror for 
HeNe; 3) Total reflection mirror for 1.054 jm (concave radius = 2.5 m); 4) Pockels cell; 5) Quarter wave plate; 6) Glan 
Taylor prism with 4 side antireflection coating; 7) Diaphragm (5.5 mm diameter); 8) Oscillator head SF 61107; 9) Output 
mirror; 10) Half wave plate; 11) Diverging lens (focal length = -600 mm); 12) Converging lens (focal length = +1000 mm); 
13) Flat total reflection mirror for 1.054 jum, 1=450; 14) Amplifier head SF 320 16; 15) Diverging lens (focal length = 
-600 mm); 16) Converging lens (focal length = +800 mm); 17) Rmax mirrors at 1.054 jm, |=450; 18) Amplifier head 

SF 320 25; 19) Electronic box of pockel cell EO 401. 


figure 9 
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Figure 9 gives a schematic overview of the laser components. The main components 
are the oscillator and two amplifiers. The oscillator consists of a cavity mirror, the am- 
plifying medium and the electro-optic trigger. The two cavity mirrors are the total re- 
flectivity concave mirror and the output mirror, which has a gaussian reflectivity profile. 
The amplifying medium is a phosphate-glass rod doped with Neodymium surrounded 
by two linear flash lamps. The two amplifiers are also equipped with a phosphate glass 
rod and have six or eight linear flash lamps, distributed around the glass rod. The two 
amplifiers receive a maximum electric energy of 8000 J each. The oscillator and am- 
plifiers are cooled by a closed circuit of deionized water. 


4.6 Data acquisition interface 


The detector outputs will be connected to adequate electronic units and, afterwards, to 
a data acquisition system that must fulfil the following requirements: (i) be fast enough 
with respect to the signal duration which is of the order of 30 nanoseconds; (ii) be 
suitable to handle the following set of measurements during each shot: (a) the detector 
outputs; (b) the incident laser power into the plasma (it is expected that the laser pulse 
power will not be entirely reproducible); and (c) the plasma background light some us 
around the laser pulse (it is advisable to make measurements before and after the 
laser pulse in order to eliminate the plasma fluctuations),. 


In the first and second phases of operation of this diagnostic (measurements at one or 
two radial positions, using 3-channels spectrometers) the detector outputs will be con- 
nected, after appropriate delays, to one or two channels of a digitalizing oscilloscope 
with the following characteristics: 250 MHz of bandwidth and 1 Gigasample of sampling 
rate for single shot operation. In the third phase of the Thomson scattering operation 
(measurements at two different radial and vertical positions, using 3-channels spec- 
trometers) and during the ISTTOK non-inductive operation (measurements with 
20-channel spectrometers), the detector output signals will be led to analog-to-digital 
conversion (ADC) modules. Each module has 12 10-bit input channels, with a mini- 
mum resolution of 0.25 pC, separate gates and conversion time of 60 us. In both 
cases, the ADC modules and the oscilloscope will be controlled by the ISTTOK control 
and data acquisition system. 


5 Conclusions 


In this paper, an outline is given of the Thomson scattering diagnostic which is going to 
be used at the tokamak ISTTOK. 


From considering the required signal-to-noise in the expected density range, we have 
concluded that a 25 J Nd:glass laser will be suitable to give information on the electron 
temperature and plasma density of the ISTTOK plasma, in a low magnetic field 


(0.45 T). Simultaneous single time measurements at two or three radial positions, us- 
ing 3-channel or 20-channel spectrometers, are envisaged. 


We hope to achieve accuracies of 1.5 to 7% for T, measurements and 10% for n, 
measurements. 


In a next paper we will present the details of the beam transport line as well as of the 
scattered light collection systems, and the calibration method for plasma density 
measurement. 
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Resumo 

Às componentes de um veclor segundo 
uma tétrada gozam da propriedade de 
serem invariantes numa transformação 
le coordenadas. Além disso, mostra-se, 
existem situações em que esses inva- 
nantes coincidem com as componentes 
tensoriais num dado reterencial de 
nércia, Nesta acepção designam-se poi 
Lsicamente mensuráveis , já que há 
rererenciais onde ganham significado 


FÍSICO. 


Virando partido dessas propriedades 
UNSCeuIu-se uma escrita universal Lud 
| as relações constitutivas entre CcOMmpo- 

es invariantes € Mostrou-se que essas 
“relações são formalmente análogas as 
que se verificam em referenciais de 


inércia, 


Abstract 


aking advantage of tetrads properties, 
entes are constructed which, besides 
being invariantom an arbitrary trans- 
tormations of coordinates also have a 
precise meaning: they coincide with 
physical quantites (felds, mductrons ) 
neasured by a given observer. In other 
words, those invariants comnede with 
neasurable quantitices. Thus, an invar 
ant and universal form of constitutivo 


lations is obtained which can be casi 


interpreted 
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CONSTITUTIVE RELATIONS I.. 
INVARIANT FORN. 


1 Introduction 


tis a consequence of the fundamental properties of tensors that they are covariant in 
a holonomic transformation of coordinates. 


In the context of Minkowski formulation, Maxwell equations are truly tensor equations 
which, in an arbitrary system of coordinates, may be written as 


DG = (1.3) 
DF" =0 (1.b) 
where D, means the covariant derivative and G'“ and “Ff are two antisymmetric sec- 


ond order tensors with six independent components each. 


In a pseudo-euclidean coordinate system each tensor component coincides with a 
physical component of a given field. Thus we have 


O H, -H, -cD, | 
-H. O H -cD, 


“E su 
| H, -H () -cD. (2.a) 
cD, cD, cD, O 
O —E, E, —cB 
Es ih E, O =, = Co 
Ee (2.b) 


-BoE O «b 
cB, cB, cB, O 


where E, H, B and D are respectively the electric and magnetic fields, magnetic and 
electric inductions measured by an observer at rest in the coordinate system. 


The components G'“ and *F* in any other coordinate system are related to Gº* and 
*FºK by means of the tensor transformation law and can be easily obtained. However 
they may not have a direct physical interpretation when considered individually or inde- 
pendently of the basis vectors and this fact should be kept in mind. 


In what concerns relations between the fields, constitutive relations, it is convenient to 
express them with the same degree of generality that it was achieved with Maxwell 
equations. 


With this purpose a universal form of constitutive relations has been searched for, that 
is, a form valid in an arbitrary system of coordinates. Once again what is at stake is to 
achieve a tensorial writing for the constitutive relations. 
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Constitutive Relations in Invariant Form 


Those authors who have a É RA tensor approach, [1], [2], begin by considering the 
tensor quantities Fu”, Gu”, u” and *Gyu* which, in fact, in the reference frame 
where the observer's Minkowski ValaBity has components u"=[0,0,0,c], coincide with E, 
cD, cB and H respectively. In that reference frame constitutive relations may be written 
as 


Gau! =ceF,u (3.a) 


Fu =cuGau (3.b) 


and, postulating that they are truly tensor equations then they will hold in an arbitrary 
coordinate system. 


This approach involves some delicate aspects. Although we shall only consider situ- 
ations where the constitutive parameters are scalars, which simpliífies the problem, 
even so some difficulties remain. In fact, constitutive relations are those traditionally 
established between physical components of the fields. In an arbitrary coordinate sys- 
tem tensorial components and physical components do not coincide nor is it possible, 
in general, to attribute a physical meaning to a tensor component, as we have already 
mentioned. Thus, a hurried interpretation of relations (3) may lead to conclusions that 
enter in conflict with usual physical notions. It is indeed sufficient to write (3) in cylindri- 
cal coordinates to verify the veracity of the previous statement. In an Euclidean coordi- 
nate system, however, everything is simple. 


For a material medium in uniform rectilinear motion it is possible to define a comoving 
reference frame where the metric is pseudo-euclidean everywhere. In that reference 
frame expressions take their simplest form and the physical meaning is immediate. 


For a medium in accelerated motion it is no longer possible to define a comoving refer- 
ence frame with pseudo-euclidean metric everywhere, although that continues to be 
possible at a single point. In an infinitesimal region around that point measurements 
cannot be distinguished from those carried out in the inertial comoving frame, although 
there exists an equivalent gravitational field. This is a fundamental step in the study of 
accelerated motion. 


In a four-dimensional space a tetrad is any set of four vectors which constitute a basis. 
Thus any tetrad may be used as a reference frame. 


It is always possible to associate a tetrad whose vectors are orthogonal and of unitary 
length with an observer in accelerated motion. On doing so an orthonormal tetrad is 
defined which provides a local, comoving reference frame with pseudo-euclidean met- 
ric, where a number of problems may be solved with elegance and simplicity, namely 
the establishment of constitutive relations. 


In section 2 the concept of orthonormal tetrad is introduced and some of the tetrad's 
properties are referred. It will be shown that the components of a tensor along a given 
tetrad do not depend on the coordinate system, and so will be called "invariants". 


In section 3 it will be shown that taking full advantage of this invariance relations can 
be found which are: 


a) universal, because valid in an arbitrary coordinate system; 
Db) consistent with (3); 
c) easily interpreted as fields measured by a given observer. 
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This subject has been included in a previous, extensive work [3]. Here we are con- 
cerned mainly with the establishment of constitutive relations which has been done in a 
straightforward manner and presenting only the results essential to our purpose. 


2 Orthonormal tetrads 


In a four-dimensional space any set of four linearly independent fourvectors is called a 
tetrad. The natural basis e, 1 to 4, defined at a given point by 


E (4) 
dx 
is a often used tetrad. 


Instead of the natural basis vectors we shall use another set of linearly independent 
| fourvectors, À(, a=1 to 4. Each one of these fourvectors is defined by its components 


along the natural basis vectors. The knowledge of the 16 components À, such as 


À 


e A, a)8i (5) 


la] 


defines the tetrad À, completely. 


Let us now consider those tetrads where each fourvector is orthogonal to the remaining 
three and of unit length. Taking the first three vectors to be spacelike and the fourth 
vector timelike, the orthonormality condition may be written as 


À i Ay ni ab) (6) 


where ".” means internal product and np; has the usual meaning nah)=0, ab, N(y= 
N(22=N(33)="1, N(44j=1. Another way of expressing the orthonormality condition is 


Saka) Kin ad (7) 
The components of a fourvector, A, along the tetrad Aa) will be called Alal(A). Thus we 
have 
A=A“(AA, = 4'e, (8) 
From (5) and (8) there results 
A'=A“(AAi, (9) 
Covariant components of À, are defined in the usual way 
À ia) — Sik Ma) (10) 
and in what concerns the indexes between parentheses the following rule is imposed 
| od E NA (11.9) 
MENA (11.b) 


where nt)=n, ap. 
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Constitutive Relations in Invariant Form 


As a consequence of (11) and the orthonormality condition (7) we obtain consecutively 


A Wu=0=0 (12.3) 
À A O, (12.b) 

and the inverse relations of (9) 
A“ (A)= ANS (13) 


The result of the contraction of two fourvectors is a scalar. Thus the components 
Ata(A) are scalars which means that they are invariant in a transformation of coordi- 
nates. 


Given two orthonormal tetrads, 4 and À, defined at the same point and related by 


ul A (14.a) 
UA es O aÃã O (14.b) 

it may be concluded from (12) that 
E jog = (15) 


This is yet another way of expressing the orthonormality condition. In fact the coeffi- 
cients that relate the Euclidean coordinates in a Lorentz transformation also verify (15). 
The indexes between parentheses are usually called "Lorentz indexes”. 


In what concerns the components of A along the tetrads u and À identical relations ap- 
pear. We have 


AS ul=E ol) (16.a) 
A ti)= La AA) (16.b) 
What has been said is easily pino to higher order tensors. We have for example 
=, ado A Pd) (17.a) | 
and its inverse 
AS, (A)= NONOAI Jo (17.b) 
as well as 


Al) = Ein, A (Ã) (17.0) 
| 


An orthonormal tetrad may be defined along the worldline of an observer in arbitrary 
motion subject to the following transport law: the timelike fourvector will be tangent to 
the observer's worldline at any given point. On doing so the components of the ob- 
servers Minkowski velocity along the spacelike vectors of the tetrad will be null and 
consequently the observer will be constantly at rest relative to the tetrad. In fact this 
problem has not one but infinite solutions, all of them related by some rotation of the 
spatial axes, [4]. 
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The tetrads just defined provide an infinite number of reference frames comoving with 
the observer where, due to the orthonormality condition, the metric Is 
pseudo-euclidean. The advantages of such reference frames have already been dis- 
cussed. 


3 Constitutive Relations 


The components of an orthonormal tetrad associated with an arbitrarily moving ob- 


server are 
=[1,0,0,0] Hilal 
=[0,1,0,0] (18.b) 
adia Hds 
=[0,0,0,1) 184) 


in the inertial frame where the observer is instantaneously at rest. 


It follows from (17.b) and (18) that the components of any tensor along À coincide with the 
tensor components along the natural basis vectors of the inertial reference frame, |. e. 


FantÃ) = Fa (19.a) 


+tubl |» mo dh 
G(A)=G (19.b) 
According to the hypothesis that local measurements cannot be distinguished from 
those made in an instantaneously comoving inertial reference frame the following 
identification is usually made, [5] 


FiatA)= Ego G=1to3 (20.a) 
FiylÃ) =cB, FostÃ) = CB, . FaghÃ) = CB; EP 
GcilAj=eDy «o t=1to3 (20.0) 


GalA)=H;, Go lA)=H,, GaylA)= H, (20.9) 


where E, H, B and D are respectively the electric and magnetic fields, magnetic and 
electric inductions measured instantaneously by the moving observer. 

If the observer is comoving with a material medium with scalar parameters £, u then 
the following relations apply 


6, o4) ( A) e CEF al À) (21 .a) 
Grag)(Ã) = fianitÃ) (21.b) 


which, due to the invariave of G (ab(A) and FrabA), are valid in an arbitrary coordinate 
system. We have thus arrived at an invariant form of the relations between fields and 
inductions. 
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Constitutive Relations in Invariant Form 


Equations (21) are consistent with (3). In fact it follows from (21) and the relations be- 
tween invariants and tensor components in an arbitrary system of coordinates, 


G; = A, A GralÃ) (22.a) 


I k 
a" FÃ) = NX Fy (22.b) 


Gu = o AO, ON tap) + cof ALA, O NaN (o = 1 Manta)E, (23) 
He É J 


Taking into account the orthonormality relations we have 


À; Na = õ! = A Na) (24) 
Substituting in (23) we obtain, after rearrangement 
l l=:“eU u 
Gu Fo Cru, + Fy u)— (25) 
Na cu c 


where u'=ch!,, is the Minkowski velocity of the medium. 


This is the form under which the constitutive relations of non-dispersive media are 
usually presented. In fact equation (25) follows immediately from (3) when the proper- 
ties of the dual tensor are considered [1]. 


Starting from (21) we have arrived at a well-known relation. Therefore it may be con- 
cluded that equations (21), derived and written in tetrad formalism are consistent with a 
tensorial form of constitutive relations and have the same degree of generality. 


We have already pointed out the difficulties in interpreting the tensor components 
along the natural basis vectors in an arbitrary system of coordinates. But the 
"invariants" that may be formed from those tensor components are easier to 
understand. However the invariants depend on the tetrad. For the observer comoving 
with the medium, observer À, the invariants are G,abj(A) and Fjapj(A) and are related by 
(21). 


Another orthonormal tetrad, 4, is associated with a second observer, in motion relative 
to the medium. The invariant components G pH) and F,apj(H) represent the fields and 


inductions measured by observer |. Since the components of a tensor along two or- 
thonormal tetrads are related by an instantaneous Lorentz transformation we obtain, 
upon substitution in (21) 


Foot du, MU) + E (uu, (1) E u | (26) 


I- ceu 


| Ê 
Gas! H ) = — Eb) | H p= 
Hc c 


where u((u) are the components of the medium Minkowski velocity. The difference 
between (25) and (26) is that in the latter expression the metric is pseudo-euclidean, 
N ab) and each invariant component has a precise meaning. 


In the final analysis equations (26) express the constitutive relations in the inertial ref- 
erence frame where observer | is seen instantaneously at rest. 


On the other hand multiplication of (25) by Ha and ao also gives (26). This is an in- 
teresting result which is, in fact, another manifestation of the invariance of the compo- 
nents of a tensor along any orthonormal tetrad. 
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4 Conclusions 


The components of a tensor along an orthonormal tetrad are scalars since they result 
from the contraction of two fourvectors. Therefore their value is independent of the co- 
ordinate system and they are abbreviately designated by "invariants”. 


This invariance under coordinate transformation would be interesting by itself and use- 
ful when applied to the study of moving media, namely electromagnetic phenomena in 
arbitrarily moving systems. As it happens the invariants along an adequately defined 
orthonormal tetrad coincide with physical quantities measured by the observer comov- 
ing with the tetrad. Also, due to the orthonormality condition, the use of tetrads allows 
an enlightening geometric interpretation. 


Invariance, physical meaning and easy geometric interpretation are three important as- 
pects, often referred in the literature, which make attractive the use of tetrads in 
Relativity. 


Based on the properties of tetrads a formalism was developed which led to an alterna- 
tive form of constitutive relations with the desirable-features of being universal and 
easily interpretable. 
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O ESCOAMENTO TURBULENTO EX 
SISTEMAS DE QUEIM- 


1 Introdução 


A estabilização de chamas de elevada intensidade em queimadores industriais e nor- 
malmente conseguida através da criação de uma zona de recirculação a jusante dos 
queimadores por indução de rotação ao ar de combustão [1]. Os processos de com- 
bustão e formação de poluentes, e em particular os fenómenos associados à estabili- 
zação e propagação das chamas, estão associados à eficiência da mistura do com- 
bustível com o ar, a qual é determinada pelas características turbulentas do escoa- 
mento na vizinhança dos queimadores. A análise de sistemas práticos de queima está 
assim dependente do conhecimento destas características, o qual é vulgarmente infe- 
rido de escoamentos sem reacção quimica tal como acontece com o desenvolvimento 
de modelos físico-matemáticos para a simulação de escoamentos com interesse em 
engenharia (e.g.[2]). As caracteristicas do escoamento turbulento em sistemas com 
combustão são no entanto alteradas pela libertação de calor dai resultante e podem, 
como consequência, ser qualitativamente diferentes daquelas tipicas de escoamentos 
não-reactivos. Exemplos incluem os processos de difusão contra-gradiente [3] em 
chamas com e sem pré-mistura dos reagentes e a turbulência gerada pelas chamas 
[4-6]. E portanto tecnicamente importante aprofundar o conhecimento dos mecanis- 
mos de produção turbulenta induzidos por combustão para optimizar o desenho e a 
operação de sistemas de queima e para validar a utilização de técnicas de simulação 
de escoamentos com flutuações de densidade (e. g. [6] e [7]), mas a realização destes 
objectivos depende da disponibilidade de dados experimentais detalhados. 


Algumas medidas experimentais das propriedades de chamas turbulentas em escoa- 
mentos recirculatórios com rotação têm sido apresentadas na literatura (e. g., [8] e 
[9)), mas a maioria destes trabalhos apenas considera valores médios de temperatura 
e velocidade e intensidades de turbulência. Por exemplo, trabalhos experimentais que 
utilizem os resultados para estimar grandezas importantes para a análise e modelação 
da turbulência, tais como balanços e correlações das tensões de Reynolds, requerem 
a determinação de gradientes espaciais de caracteristicas vectoriais e não têm sido 
apresentados com o detalhe suficiente para possibilitar a compreensão dos mecanis- 
mos de mistura e transporte turbulento em escoamentos com interesse em engenharia 
(e. 9., [10]). Só muito recentemente a velocimetria laser tem sido utilizada para obter 
este tipo de informação (e. g., [5] e [7]), embora em escoamentos relativamente sim- 
ples, tais como camadas de mistura sem reacção quimica, jactos axissimétricos com 
combustão e com niveis moderados de rotação, ou chamas estabilizadas em corpos 
não fuselados. Os escoamentos recirculatórios com elevados níveis de rotação típicos 
dos sistemas de queima industriais estao associados a processos resultantes da acção 
de forças centrifugas e do efeito da curvatura das linhas de correntes no escoamento 
turbulento em regiões de elevadas tensões de corte (e. g., [6] e [1]), e não têm sido 
suficientemente documentados. O objectivo deste trabalho consiste portanto em con- 
tribuir para a compreensão dos processos de mistura turbulenta em escoamentos tipi- 
cos de sistemas de queima, a qual é essencial para se aperfeiçoar o projecto de dis- 
positivos práticos. Sao apresentadas medidas detalhadas dos campos de velocidade, 
de temperatura e concentrações de especies quimicas do escoamento na vizinhança 
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de um modelo laboratorial de um queimador industrial e são estimados os termos mais | 
importantes que contribuem para a conservação de energia cinética turbulenta. O mo- 

delo experimental é descrito na secção seguinte e os resultados são apresentados e 
discutidos nas Secções 3 e 4, respectivamente. A ultima secção sumariza as princi- 
pais conclusões deste trabalho. 


2 Configuração do Escoamento e Método Experimental 


A instalação experimental utilizada, figura 1, é essencialmente constituida por um 
atomizador comercial (Delavan SNA 1.00) com 23 mm de diâmetro rodeado por um 
escoamento coaxial de baixa velocidade de combustivel gasoso (diâmetro interi- 
or=54mm). o qual está por sua vez colocado num escoamento coaxial exterior de alta 
velocidade (diâmetro interior=84 mm). A instalação inclui ainda uma secção divergen- 
te ("quar!' na literatura inglesa), tipica das que se encontram em queimadores indus- 
triais, e que pode ser removida para a medição das condições de fronteira. Os escoa- 
mentos interior e exterior são animados de rotação através de dois conjuntos de pás 
planas fixas, inclinadas a 45 . resultando em numeros de rotação respectivamente 
iguais a 5,=0.77 e 50.85 nos escoamentos exterior e interior. O caudal de ar é me- 
dido com um orifício calibrado com uma precisão de 2%, enquanto o caudal mássico 
de propano é medido por um rotâmetro com um erro absoluto inferior a 0.1 g/s, cor- 
respondendo a uma precisão de 1.75 % para as condições utilizadas neste trabalho. 


| 117 mm | 


40 mm 


| ee 54 mm | 


Do = 84 mm 


figura 1 Diagrama do queimador laboratorial 


As medidas aqui apresentadas foram obtidas nos escoamentos com e sem combus- 
tão, sem injecção de combustive! liquido através do atomizador. No escoamento com 
combustão as velocidades médias, definidas como o quociente entre o caudal volúmi- 
co e a area da secção transversa|, correspondem a U,=30 m/s (Re,=49500) no esco- 
amento exterior e U,;s=1.8 ms (hi =3000) no escoamento interior, o que corresponde 
a uma chama com uma razão ar/combustível igual a 27.6 e uma taxa de libertação de 
energia igual a 350 kW. O escoamento sem combustão corresponde às mesmas con- 
dições de operação, tendo o combustivel gasoso sido substituído por ar com o mesmo 
fluxo de quantidade de movimento, correspondendo a um numero de Reynolds igual: a 
4000 no escoamento interior. 
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ESCOAMENTO SEM COMBUSTÃO 
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figura 2 
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ESCOAMENTO COM COMBUSTÃO 


Distribuição espacial de linhas de corrente para os escoamentos com e sem combustão 
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O queimador está montado na posição vertical numa mesa de atravessamento com 


plano de saida do queimador e na linha central do mesmo. 
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figura 3 
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deslocamento nas três direcções ortogonais, permitindo o seu posicionamento com 
uma precisão de +0.25 mm. A origem do sistema de coordenadas (x,r,8) foi tomada no 


Caracteristicas aerodinamicas dos escoomentos médios com e sem combustão: linha de corrente nula («p = ()); linha de 


velocidade axial nulo (! =) linha de velocidade axial maxima (6)! tor = (1) fronteira do escoamento ((). E. ). 
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As medidas de velocidade foram efectuadas com um velocimetro laser de uma com- 
ponente o qual inclui uma fonte de luz laser de Argon-lon com uma potência nominal 
de 1 W a funcionar no comprimento de onda da luz verde (514.5 nm) e sensibilidade à 


direcção do escoamento através de modulação óptico-acústica (células de Bragg). O 
sistema de transmissão de luz incluía alternativamente uma lente de 310 mm ou de 
600 mm de distância focal de modo a manter o sistema óptico afastado da chama, re- 
sultando em volumes de medição com eixos ortogonais iguais a 0.135 mm e 
1.527 mm ou 0.262 mm e 5.074 mm respectivamente. A luz dispersa por partículas de 
alumina adicionadas aos escoamentos de ar e combustível era colectada em “forward 
scatter" e detectada por um fotomultiplicador. Os sinais Doppler daí resultantes foram 
processados por um contador de frequência comercial e são considerados correspon- 
derem a valores ponderados com a densidade (e.g., [12]). Os detalhes do velocímetro 
e da precisão das medidas foram no entanto descritos por Moreira [13] e não são aqui 
discutidos. 


Os valores médios de temperatura foram obtidos com termopares de Pt'Pt13%Rh fta- 
bricados com fios de 0.040 mm de diâmetro. Os fios eram suportados por uma sonda 
construída com o mesmo material mas com fios de 0.500 mm de diâmetro, a qual foi 
introduzida no escoamento de modo a minimizar os efeitos associados com a intrusão 
da sonda. Uma sonda de aspiração de aço arrefecida a água com um diâmetro exteri- 
or de 7 mm foi utilizada para extrair amostras locais dos gazes de combustão de modo 
a medir valores médios ponderados com a densidade das concentrações de COs, CO, 
O,, e hidrocarbonetos não queimados, UHC. Os resultados foram utilizados para de- 
terminar a distribuição local de fracção de mistura, f(i.e., fracção mássica de combus- 
tível queimado e não-queimado), a qual é analisada nas secções seguintes juntamente 
com os campos de velocidades e temperatura. 
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figura 4 Distribuição axial do caudal máximo recirculado nos escoamentos com e sem combustão. 
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3 Resultados 


A figura 2 apresenta a distribuição de linhas de corrente dos escoamentos com e sem 
combustão e permite observar as características mais relevantes do escoamento mé- 
dio na vizinhança da saída do queimador. Os dois escoamentos apresentam caracte- 
rísticas semelhantes e típicas de escoamentos com elevados níveis de rotação, e.g. 
[14-16], que incluem uma zona de recirculação central induzida pela rotação, rodeada 
de uma região de forte aceleração axial com elevados gradientes radiais e onde a 
componente tangencial da velocidade local atinge o seu valor máximo. Apesar das 
semelhanças qualitativas entre os dois escoamentos, a combustão induz diferenças 
quantitativas significativas, como apresentado nas figuras 3 e 4: as componentes axial 
e tangencial da velocidade média aumentam devido à diminuição da densidade, uma 
vez que a quantidade de movimento se mantém constante; o caudal mássico recircu- 
lado diminui de cerca de 67% para 14% do caudal de ar insuflado; o comprimento da 
zona de recirculação diminui cerca de 32.5% e a sua largura máxima aumenta cerca 
de 12%. 


FRACÇÃO DE MISTURA, f 


r (mm) 


TEMPERATURA, T 


figura 5 Distribuição espacial de linhas de valor constante de fraccão de mistura, f , e de temperatura média, 7. 


Os resultados da figura 5 mostram que a zona do recirculação do escoamento com 
combustão corresponde a uma região de elevadas temperaturas à semelhança do que 
acontece em outras chamas estabilizadas em zonas de recirculação (e.g., [7,8]). A lIi- 


nha de estequiometria local (i.e., / = 0.060) é exterior à linha de velocidade axial nula 
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(ver figura 3) e, consequentemente, a zona de recirculação está associada à existên- 
cia de excesso de combustível. As concentrações de monóxido de carbono e hidro- 
carbonetos não queimados anulam-se no interior da zona de recirculação, enquanto 
que ao longo da camada de corte os níveis de CO atingem 3.5%. Isto resulta do ba- 
lanço entre a formação inicial de CO devido à oxidação de CaHg e à oxidação subse- 
quente de CO em CO, no interior da zona de recirculação. 


A figura 6 apresenta as distribuições de energia cinética turbulenta, k, e de tensões de 


do 
Reynolds uv” dos escoamentos com e sem combustão e mostra que a reacção qui- 
mica não altera significativamente o campo turbulento na região a montante de 
x/Dg=1.0, onde as flutuações de velocidade têm valores pequenos no interior da zona 
de recirculação e elevados na camada de corte anelar onde as taxas de deformação 
são elevadas. A jusante de x/Do=1.5 ambos os escoamentos apresentam característi- 
cas distintas: na ausência de reacção química a energia cinética turbulenta apresenta 
valores máximos ao longo da camada de corte adjacente à zona de recirculação e as 
tensões normais diminuem com a distância à saída do queimador; por outro lado, no 
escoamento com combustão os valores máximos de k ocorrem na vizinhança do 
ponto de estagnação, onde v”” e w”” são consideravelmente superiores a u”, à se- 


melhança do que foi observado, por exemplo, por [7] e [17] em escoamentos com re- 
circulação introduzida por corpos não fuselados. 
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figura 6 Distribuição espacial de linhas de valor constante de energia cinética turbulenta, k, e tensão de Reynolds uv" para os 
escoumentos com e sem combustão. 
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Os níveis de anisotropia tipicamente encontrados para as tensões de Reynolds no es- 
coamento com combustão podem ser facilmente identificados na figura 7, na qual são 
apresentados perfis radiais das três tensões normais na zona de máxima largura da 
zona de recirculação e perto da zona de estagnação. Nesta última região, para 


x/D,=2.5, os valores consideravelmente elevados das flutuações turbulentas das 


componentes transversais de velocidade estão associados a flutuações periódicas do 
escoamento com uma frequência de 30 Hz. O número de Strouhal correspondente, 
com base na definição de [14] e [18] para escoamentos com rotação, é igual a 0.2 e, 
consequentemente, uma ordem de grandeza inferior aos valores tipicamente encon- 
trados na literatura para caracterizar oscilações de zonas de recirculação em torno do 
seu eixo de simetria (i.e., “precessing vortex core” na literatura inglesa). Por outro 
lado, se a frequência de oscilação observada no presente caso for adimensionalizada 
com base nas escalas de velocidade e comprimento características da camada de 
corte que envolve a zona de recirculação, o número de Strouhal resultante é típico dos 
processos de "bursting" associados a oscilações de grande escala em camadas de 
corte, como descrito por [19] e [26] e analisado em detalhe por [13]. O fenómeno é no 
entanto uma característica da presente geometria e a sua extrapolação para outros 
sistemas de queima não deve ser efectuada directamente. Deve ser ainda notado que 
a energia associada às flutuações periódicas é apenas 4% da energia total contida no 
espectro de flutuações, pelo que a camada de corte envolvente da zona de recircula- 
ção deve ser considerada como uma região de intensa geração de turbulência. 


a) 30.0 pr 


20.0 


figura 7a Perfis radiais das tensões normais para o escoamento com combustão. 


Perfil radial na zona de máxima largura da bolha de recirculação a x /D; = 1.0. 


figura 7b Perfis radiais das tensões normais para o escoamento com combustão. 
Perfil radial na zona de estagnação, x/D; = 2.1. 


à ie 
As distribuições de tensões de Reynolds wu v (figura 6) apresentam valores máximos 


4 


coincidentes com os de u””, com valores mais elevados no escoamento com com- 
bustão. Para os escoamentos aqui considerados o sinal da tensão de corte está relaci- 
onado com o da taxa de deformação dU /dr de acordo com a hipótese de viscosidade 
turbulenta (e. g., [17]), excepto na região inicial da camada de corte adjacente à zona 
de recirculação no escoamento com combustão, onde a taxa de deformação é pró- 
xima de zero. Os valores máximos do coeficiente de correlação da tensão de 
Reynolds aumentam cerca de -0.3 à saída do "quarl" (i.e., x/Dqg = 0.6) até cerca de - 
0.7 em x/Do = 1.0 e, posteriormente, diminuem para valores constantes e iguais a - 
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0.4. Resultados semelhantes foram apresentados por [15] na zona de recirculação de 
uma chama de pré-mistura estabilizada num disco, tendo sido atribuídos à presença 
de taxas adicionais de deformação associadas à curvatura das linhas de corrente (ver 
também [24]). 


a) 1.0 b) 10 
0.5 
s 
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É 
- 
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figura 8 Perfis radiais das velocidades de transporte turbulento no escoamento com combustão. 


Os momentos de ordem superior das flutuações de velocidade estão associados aos 
processos de transporte turbulento e têm sido analisadas com base na acção convecti- 
va das grandes escalas turbulentas, e. g., [22] e [23]. No presente trabalho foram 


! sra 


efectuadas medidas dos momentos de terceira ordem uv", v” e uv” que repre- 


: Pro Pr 


sentam, respectivamente, os fluxos radiais de u”, v” e w'v”. Os resultados obtidos 
para qualquer dos escoamentos estudados mostram que o transporte turbulento das 
tensões normais e de corte se dá no sentido dos respectivos gradientes radiais e apre- 
sentam desenvolvimentos semelhantes aos obtidos por outros autores a jusante de zo- | 
nas de recirculação, e. g. [24], [25] e [26]. A figura 8 mostra a componente radial das 


velocidades de transporte da tensão de corte w'v” e da energia cinética turbulenta, k, 
definidas como em [27] por: 


pu ET | 
FUV a | 
VP MV | 


e sugere que na camada de corte adjacente à zona de recirculação (onde JU /dr >0) 
os valores não são afectados pela presença de combustão: as velocidades de trans- 


porte são negativas na região interior e positivas na região exterior, sendo nulas nos 
picos (negativos) da tensão de corte e da energia cinética turbulenta. Na camada de 
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corte exterior a razão Vyyw'Vk é cerca de 2 e, portanto, inferior ao valor de 3.6 caracte- 
rístico de camadas de mistura planas ([28] e [29]), mas próximo do valor de 2.5 obtido 
por [21] numa camada limite curva. A diminuição de V,, e V, para jusante é consisten- 


te com a diminuição da actividade turbulenta. No entanto, a tensão de corte uv” na 
camada exterior depende mais do transporte turbulento do que a manutenção de k 
(uma vez que V,,> Vi), significando que a tensão de corte é mais sensível a uma redu- 
ção da velocidade de transporte do que a energia cinética turbulenta. 


4 Discussão 


Os resultados analisados nos parágrafos anteriores quantificam o efeito da combustão 
nas características aerodinâmicas de um sistema de queima típico e identificam as- 
pectos importantes da termoquimica destes sistemas com aplicações directas na sua 
optimização no que respeita à eficiência de combustão e à minimização dos poluentes 
gerados pela combustão. Verifica-se que no interior da zona de recirculação vulgar- 
mente associada a queimadores industriais o combustivel é consumido e o CO oxida- 
do em COs. A zona de recirculação é essencialmente constituída por produtos da 
combustão cuja composição não satisfez as condições de equilibrio quimico. Os valo- 
res máximos de CO e de hidrocarbonetos não queimados ocorrem na camada de 
corte em zonas de estequiometria local e, portanto, os resultados sugerem que as limi- 
tações de ordem química impostas na queima do combustível e de CO devem ser 
desprezáveis. No entanto, a estabilização de chamas em sistemas práticos de com- 
bustão depende da distribuição de estequiometria local (e.g., [1]), a qual é função dos 
processos de mistura turbulenta entre o combustivel e o oxidante e, portanto, é influ- 
enciada pelas caracteristicas turbulentas do escoamento. Os resultados aqui apresen- 
tados incluem informação detalhada sobre estas caracteristicas e possibilitam avaliar 
a precisão de métodos de cálculo do escoamento em sistemas de queima de modo a 
interpolar e extrapolar a informação adquirida experimentalmente. 


Os resultados mostram que no escoamento aqui em estudo os níveis de turbulência 
são amplificados pela combustão, apesar das velocidades médias também aumenta- 
rem em resultado da reacção quimica sem que se observe um correspondente au- 
mento significativo da espessura das camadas de corte. Consequentemente, é de es- 
perar um aumento dos processos de geração de energia cinética turbulenta por inter- 
acção das deformações com as tensões de corte (e.g., [30]). A questão que se coloca 
está associada à existência de mecanismos de geração de turbulência para além da- 
quele resultante do aumento da taxa de deformação, tais como sejam a interacção 
entre as flutuações de densidade e a presença de gradientes de pressão média, e. g., 
[6]. [7] e [111]. 


As medidas de velocidades média e turbulenta apresentadas anteriormente permitem 
estimar os termos de convecção e produção nas equações de transporte de quantida- 
de de movimento e da energia cinética turbulenta, os quais são aqui usados para ana- 
lisar os mecanismos envolvidos na geração de turbulência em sistemas típicos de 
queima de combustíveis. As estimativas são aproximações devido ao erro associado à 
determinação dos gradientes espaciais, mas os valores obtidos são suficientemente 
precisos para possibilitarem analisar a importância relativa dos vários termos envolvi- 
dos nas equações. Para escoamentos com elevados números de Reynolds a equação 
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de transporte das tensões de Reynolds na forma de valores médios ponderados com a 
densidade é dada por: 


“ auu] E :9U;| 
PU, o — pu u; Ou; Ah 
“ax, dx, ox, ] 

ó pu, dd: My — HU; - OP +, a 
“9x Ox, “9x, 
=p = 
Au, Eu e 
“A 'j dx 


(14) 

A figura 9 mostra perfis radiais dos termos convectivos, |, e de produção, Il, da equa- 
ção (3) obtidos na zona de máxima largura da bolha de recirculação e perto da zona 
de estagnação dos escoamentos com e sem combustão, estando os resultados repre- 
sentados de modo a que valores negativos de convecção representam um ganho de 
energia cinética turbulenta. Os termos de produção turbulenta incluem a interacção 
quer de tensões normais, quer de tensões de corte, com taxas de deformação, es- 
tando os dois processos devidamente diferenciados na figura. Os termos de difusão, 
Ill, e dissipação turbulenta não estão representados por limitações experimentais e, 
juntamente com o termo de produção associado à acção dos gradientes de pressão 
média, IV, foram estimados de modo a satisfazer o balanço dos outros termos. 


a) Escoamento Isotérmico 
100.0 


x/D, =2.14 


bi 


30.0 DEDE 10.0 
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figura 9 — Pertis radiais dos termos da equação de transporte de energia cinética turbulenta nos escoamentos com e sem combustão. 
O - Produção por tensões normais 4 - Produção por tensões de corte CL - Convecção — - Diferença 
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Os resultados mostram que a zona de recirculação corresponde a uma região de fraca 


actividade turbulenta e que na camada de corte exterior (aonde OU /dr<0) a distribu- 


ição dos vários termos é semelhante à encontrada em camadas de mistura não afec- 
tadas por efeitos de curvatura das linhas de corrente (e.g., [24] e [31]): A produção de 
k é essencialmente devida à interacção entre tensões e deformações de corte e é ba- 
lanceada por processos de difusão e dissipação turbulenta. Por outro lado, no centro 
do jacto anelar a produção turbulenta por tensões normais e de corte é desprezavel e 
os termos convectivos são os maiores representando um ganho de k. Interiormente ao 
jacto anelar e na proximidade da linha de corrente de separação, a distribuição dos 
vários termos para ambos os escoamentos estudados assemelham-se aos resultados 
obtidos por [24] numa camada de corte de mistura e por [24] e [26] a montante de um 
zona de recolagem após separação do escoamento: os termos convectivos são pe- 
quenos, a produção turbulenta é essencialmente por tensões de corte e é balanceada 
pela difusão e dissipações turbulentas. Perto de x/D,=2.1, a deflexão das linhas de 
corrente para os dois escoamentos estudados (ver figura 1) ocorre numa zona onde a 
produção turbulenta pela interacção entre tensões e deformações normais é importan- 


Al na 
te. Este facto é consistente com o aumento da razão v” / u”” observada na zona de 


estagnação e está de acordo com as observações de [7] numa chama turbulenta es- 
tabilizada num disco. 


As distribuições da figura 9b) não identificam a contribuição dos termos fonte devidos 
aos gradientes de pressão média na conservação de energia cinética turbulenta, os 
quais estão associados a geração de turbulência pela chama como discutido anterior- 
mente. É evidente que a importância relativa destes termos depende da coexistência 


se 


de valores elevados de gradientes de pressão dp/dr e de w, , os quais dependem das 


condições particulares em estudo. Por exemplo, [6] mostraram que estes termos são 
desprezáveis em chamas sem pré-mistura tipo jacto, mas [32] consideram que os 
mesmos termos podem ser superiores aqueles que envolvem a interacção de tensões 
e deformações de corte numa chama confinada com rotação. Adicionalmente, [7] 
mostraram que numa chama com pré-mistura de reagentes estabilizada num disco 
estes termos são desprezáveis na conservação de k, mas constituem fontes importan- 
tes na conservação dos fluxos turbulentos de calor dando origem a processos de difu- 
são “contra-gradiente". A análise dos termos das equações de conservação de mo- 
mento permite estimar os valores dos gradientes de pressão, e para a chama aqui 
estudada (ver detalhes por [13]), dp/dx e dp/dr são da mesma ordem de grandeza na 
vizinhança da zona de estagnação, mas têm sinais contrários e actuam respectiva- 


mente como poço de u” e fonte de vo . Apesar de não ser possivel concluir que os 
dois termos se cancelam mutuamente como em [33], é evidente que contribuem para 
o aumento da anisotropia das flutuações de velocidade observadas anteriormente. A 
implicação directa destes resultados para a modelação fisico-matemática do escoa- 
mento em sistemas de queima significa que é necessário considerar equações de 
transporte individualmente para cada tensão normal e que os processos relacionados 
com a ocorrência de flutuações de densidade têm de ser representados directamente 
se se pretender simular correctamente os detalhes da aerotermodinâmica dos escoa- 
mentos. 
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5 Conclusões 


Medidas de velocidade, temperatura e concentrações de espécies químicas efectua- 
das na vizinhança de um modelo laboratorial de um queimador industrial, permitiram 
estudar os efeitos da combustão nos processos de transporte turbulento típicos dos 
observados em chamas estabilizadas em zonas de recirculação induzidas por rotação. 
O estudo envolveu a utilização de um velocímetro laser, termopares finos de Pt/Rh e 
de sondas de aspiração juntamente com um sistema dedicado de análise de gases da 
combustão, e refere-se à queima de propano gasoso para uma potência térmica de 
350 KW. Os resultados mostram que o escoamento na vizinhança do queimador é 
dominado pela presença de uma zona de recirculação de grandes dimensões que cor- 
responde a uma região de produtos da combustão a altas temperaturas. Os campos 
de velocidade dos escoamentos com e sem combustão são qualitativamente seme- 
lhantes, mas a combustão diminui o caudal mássico recirculado e o comprimento da 
zona de recirculação e aumenta a sua largura máxima e, portanto, a curvatura das li- 
nhas de corrente. Estas características impõem efeitos importantes no campo turbu- 
lento, o qual é anisotrópico apresentando desvios acentuados de isotropia na proximi- 
dade da zona de estagnação. As distribuições temporais de velocidade nesta zona 
evidenciaram a presença de flutuações de grande escala com uma frequência de 30 
Hz, apesar da sua contribuição para o valor global da energia cinética turbulenta ser 
inferior a 4%. 

A analise dos termos das equações de transporte das tensões de Reynolds mostra que 
a interacção entre as tensões normais e as taxas de deformação normal influenciam 
fortemente o campo turbulento na vizinhança do ponto de estagnação e sugerem que 
taxas de deformação adicionais, tais como aquelas associadas ao campo de pressão 
médio, promovem pelo menos a anisotropia das tensões de Reynolds. Os resultados 
sugerem ainda que a difusão e dissipação de energia cinética turbulenta são compara- 
tivamente importantes, em particular na região de estagnação e no jacto anelar com 
rotação. 
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O caminho para a protecção ambiental integrada 


HM A Hoechst tem sido pioneira no tratamento de efluentes líquidos. Em meados do anos 60: 
instalámos a primeira estação de tratamento biológico de esgotos industriais da Europa. 
Agora, já na terceira geração, o ºBiohoch-Reaktor desenvolvido em cooperação com a 
nossa associada UHDE representa a última palavra da tecnologia a nível internacional. 

aih Is == 7 E da purificação de efluentes gasosos. 
7 ao ir D/A Ro É tn A protecção ambiental hoje integrada no 
DL di O TE E4é centro das atenções da nossa investigação e 
| ns à te dE desenvolvimento, tem em vista processos de 

produção que poupem matérias primas e 

energia, e minimizem ou, no caso ideal, 

anulem as emissões poluentes, bem como 
processos para a recuperação de matérias 
primas a partir de resíduos. 


baterias herméticas, separador de fibra de vidro, 
bateria de chumbo-ácido 


Resumo 

O separador de microfibra de vidro é o 
nico empregue nas baterias herméticas 
le recombinação de gás (5LA), Apesar 
de ser dotado de boas propriedades 
químicas e eléctricas, possui fraca 
onsistência mecânica e é um poluente, 
elo que se procura obter alternativas a 
este separador. Neste artigo discute-se 
sta problemática, propondo-se alterma- 


tivas é analisando estudos em curso, 


Abstract 

we problematic of hattery separators for 
SLA batteries 15 examined. A viable 
altemative to existing polymene enve- 
lope type materials 15 considered. 


sealed batteries, battery separators, lead-acid battery 
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PROBLEMÁTICA DO SEPARADOR 
NUMA BATERIA HERMÉTICA DE 
RECOMBINAÇÃO DE GAS (SLA) 


| Generalidades 


No início dos anos 70 surgiu a nova tecnologia de bateria chumbo-ácida selada, re- 
combinante de gás (SLA.) criada à semelhança do tipo níquel-cádmio hermético, fun- 
damentada num mecanismo interno de recombinação de oxigénio nas placas negati- 
vas. Esta recombinação produz água, reduzindo drasticamente o seu consumo e pos- 
sibilitando a realização duma verdadeira bateria sem manutenção. 


Desde essa data, assiste-se a uma procura sempre crescente (4-5% anual) deste tipo 
de bateria que tende a substituir as tradicionais baterias de placas Planté, tubular e/ou 
plana em aplicações estacionárias (telefones, telecomunicações, luz de emergência, 
centrais de potência, alimentação ininterrupta de sistemas computorizados especial- 
mente da gama de 20 a 100 ah/10). A bateria SLA está projectada para uma vida 
média de 10 anos. Espera-lhe um futuro brilhante, consumo em massa em bateria de 
automóvel, se forem ultrapassados problemas de fiabilidade e custo. 


2 Aspectos Construtivos duma Bateria SLA. 


Dum modo resumido, podemos referir: 


2.1. Armaduras negativas à base de chumbo puro ou mais frequentemente à base de 
liga de chumbo-cálcio-estanho, com grande sobretensão para a libertação de hi- 
drogénio (dificulta a libertação de hidrogénio). Dimensionadas para comportar ex- 
cesso de matéria activa, a fim de em operação à carga, não ocorrer libertação de 
hidrogénio, mas simplesmente oxigênio para o qual, o elemento dispõe de meca- 
nismo de recombinação. 


2.2. Armaduras positivas, geralmente isentas de antimóônio (também em chumbo puro 
ou liga chumbo-cálcio-estanho) para evitar a contaminação da pasta negativa 
(com antimónio) o que iria facilitar a libertação de hidrogénio. 


2.3. À escolha do separador é vital para a eficiência de recombinação do oxigênio no 
interior da bateria. O separador consiste geralmente num feltro de microfibra de 
vidro, feito pelo processo tradicional de fabrico de papel, extremamente puro e 
com elevada porosidade; nele está embebido todo o electrólito, numa condição 
não completamente saturada de modo a criar canais de gás para a passagem do 
oxigénio (libertado à carga na placa positiva) para a placa negativa. 


2.4. É fundamental assegurar a estanquicidade permanente da bateria, pois o acesso 
de ar ambiente descarregaria completamente a placa negativa e ocorria perda de 
água do sistema por evaporação. 


2.5. À bateria tem uma válvula de segurança, de sentido único, que abre quando ocor- 
re sobrepressão de 200-350 gem"? (sobrecargas ocasionais). 


Original submetido para publicação em 10 de Setembro de 1991 
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2.6. O monobloco e tampa são moldados em plástico (geralmente ABS.), com módulo 
de elasticidade e resistência à fractura elevados, de modo a não deformarem em 
serviço. Também possuem resistência a propagação de chama. 


3 Mecanismo de Recombinação 
As reacções electroquimicas de carga-descarga são as mesmas nos 2 tipos de bateri- 
as: 
Herméticas, recombinantes. 
- Convencionais, de electrólito livre. 
Na placa positiva, as reacções podem ser expressas por: 
PbO,+ H,S0,+2H' +2€ & PbS0,+2H,0 

Na placa negativa são expressas por: 

Pb+SO <> PBSO, + 2e 
Globalmente, teremos: 

PbO,+ PDA ZH,5SO, 6 2PbSO,+2H,0 


Na parte final da carga, juntamente com a reconversão das matérias activas nas pla- 
cas, ocorre decomposição de água em oxigénio e hidrogénio que escapam para a at- 


mosfera, | 
H,0>H+- 


2] 

O, 

Neste aspecto, a bateria hermética estrutura-se, para não libertar hidrogénio e recom- 
binar oxigénio. Evita libertar hidrogénio na placa negativa sobredimensionando esta 
em relação à positiva. 


A placa negativa encontra-se sempre num estado parcial de descarga embora muito 
próximo de 100% carga. O hidrogénio só começa a libertar-se, quando a placa negati- 
va atinge 98%. 


Recombina o oxigénio segundo a seguinte sequência de reacções: 
| 
a) Pb+> O = PhO 
b) 2PbO+2H,SO, — 2PbSO,+2H,0 


co) PbSO,+2H'+2€ > Pb+HSSO, 


Para isso todo o oxigénio gerado na positiva, não deve escapar para a parte superior 
acima das placas (na bateria hermética ao contrário da convencional, não tem câmara 
de gases nem câmara de sedimentos). 


O oxigénio deve difundir-se rapidamente através de canais gasosos, estruturados no 
separador de microfibra de vidro, da positiva para a superfície da placa negativa. 


Esta condição de canais gasosos através do separador, é criada impondo saturação 
parcial deste, com electrólito (e não saturação completa como ocorre na bateria con- 
vencional). A difusão do oxigénio no electrólito é muito lenta, devido à diminuta absor- 
ção daquele gás no electrólito. Pelo contrário, a difusão do oxigénio no ar, é relativa- 
mente rápida. 
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A extensão da reacção a) é condicionada pela velocidade de difusão do oxigénio 
nos canais gasosos do separador e é esta que condiciona a eficiência de 
recombinação. 


Em termos práticos, esta limitação condiciona o valor máximo de intensidade de cor- 
rente de carga ou de tensão a carga. 


Se ultrapassarmos estes valores, parte do oxigénio não tem possibilidade de reagir 
com a matéria activa das placas negativas (Pb esponjoso) e liberta-se para a parte 
superior acima das placas. À pressão interior aumenta acima da regulada na valvula 
de segurança; esta abre e descarrega o gás para a atmosfera, ocasionando uma perda 
significativa de eficiência. A bateria “seca um pouco” e perde capacidade eléctrica. 


Em termos gerais, os fabricantes de baterias limitam a corrente máxima de sobrecar- 
ga a C,/20A (em que C,, é a capacidade nominal no regime de 10 horas) e/ou a ten- 
são máxima de carga a 2,5 V/Elemento. 


Isto não tem em conta: 


a) Não produzir oxigénio a uma taxa superior à permitida pela velocidade de acesso 
daquele gás à superficie da placa negativa. 


b) Limitar a taxa de libertação de calor (reacção de oxidação 420, + Pb — PbO) de 
modo que a bateria não atinja valores de temperatura elevada em serviço. A bate- 
ria tem sempre dificuldade em escoar o calor gerado no seu interior. 


Como referimos atrás, o separador desempenha papel primordial neste transporte ga- 
soso. É desejável escolher o material e estrutura de poros que favoreçam no máximo 
a velocidade de difusão do oxigénio. A microfibra de vidro apresenta propriedades 
notáveis para esta aplicação e é usada extensivamente só ou misturada com outras 
fibras plásticas hidrofiílicas (poliéster por exemplo). 


Note-se que os compostos de sílica , como se sabe, desempenham uma importância 
primordial no transporte de oxigênio da placa positiva para a placa negativa. 


4 Desenvolvimento de Novos Separadores Poliméricos 


O separador de microfibra de vidro é o único empregue nas baterias seladas que es- 
tão presentemente a ser comercializadas; isto deve-se a algumas das suas proprieda- 
des (resistência química a oxidantes enérgicos - peróxido de chumbo e oxigénio livre 
em ácido sulfúrico diluído, d=1,300; baixa resistência eléctrica; elevada capacidade de 
absorção do electrólito). Este separador é obtido pelo processo de via húmida usado 
na industria do papel. Apesar de ser dotado de boas propriedades, obriga a cuidados 
especiais aquando da sua montagem nas baterias seladas tipo recombinante, dada a 
sua fraca consistência mecânica. Tem também outros inconvenientes, nomeadamente 
o de ser um poluente, pelo que se procura obter alternativas a este separador. 


Nessa medida, iniciou-se uma pesquisa bibliográfica, adquiriram-se livros, artigos e 
patentes e, da análise deste material, concluiu-se que o separador alternativo (maior 


flexibilidade, resistência à abrasão, ...) poderá ser um feltro, não tecido, composto por 
lá de vidro e fibras sintéticas (polipropileno, polietileno, poliéster, poliestireno, poliacri- 
lonitrilo) extremamente finas (1 a 10 |), obtido por tratamento térmico (graft) e mecã- 
nico (compressão). Nesta estrutura o polipropileno actuaria como reforço mecânico e 
o custo de fabricação do separador seria diminuído (na produção de fibra de vidro há 
elevado consumo de energia). Posteriormente haverá que sujeitar o feltro a um trata- 
mento de humedecimento pelo electrólito (solução de ácido sulfúrico). Esta estrutura 
deverá ter cerca de 85 - 95% de porosidade em volume e um tamanho de poro máxi- 
mo inferior a 30 n. 
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Uma vez estabelecida uma primeira ideia acerca do separador que propomos obter e 
vir a utilizar, foram contactadas várias firmas nacionais (Quimigal, Fisipe, Finicisa, 
Ipetex, Geodreno Geotexteis, Fanafel, Tent e Portucel) e estrangeiras (American Felt 
& Filter Company, Johnson-Progress Ltd. , Texon, Hércules) fornecedores de filtros e 
feltros. Estas firmas têm-nos cedido amostras, que têm sido utilizadas em filtros de 
prensa e com outros fins, e cuja caracterização está em curso. 


Simultaneamente iniciaram-se ensaios preliminares de elementos de acumuladores 
adquiridos às firmas Oldham, Hitachi, Absolyte e Tudor Espanhola. 


Os resultados destes ensaios e doutros estudos prosseguidos neste dominio farão 
parte dum trabalho a publicar oportunamente. 
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Normas de Publicação de Artigos na Técnica, 
revista de Engenharia 


1 A revista TÉCNICA, em princípio, publicará 
anualmente quatro números de natureza técni- 
co-científica. 


2 Atendendo ao contexto onde se insere e público 
a que, em conformidade e de modo principal, se 
destina, a revista integra-se a partir das três se- 
guintes formas de publicação: 


ARTIGOS POR CONVITE 


São da iniciativa da Direcção da Técnica e desti- 
nam-se, em princípio, a trabalhos de síntese no 
domínio das ciências básicas ou ciências da en- 
genharia. Além do mais, procura-se que estes ar- 
tigos se constituam fontes recomendadas de es- 
tudo aos alunos e de consulta para investigado- 
res e docentes. 


ARTIGOS TÉCNICO-CIENTÍFICOS 


São artigos submetidos a publicação. Desti- 
nam-se essencialmente a dar noticia sobre aqui- 
sições actuais da ciência e da técnica mas em 
termos que se vá além da simples divulgação, 
sendo desejável que nesses trabalhos figurem 
parâmetros e referências que permitam elaborar 
juízos críticos e comparativos. Fundamentalmen- 
te, uma notícia sucinta, com carácter crítico, dan- 
do conta de uma aquisição recente e importante 
para a técnica é a marca fundamental destes arti- 
gos que vão traduzir, em princípio, investigação 
científica realizada pelo autor. 


NOTAS CIENTÍFICAS 


Destinam-se essencialmente à apresentação de 
novos resultados no domínio da investigação ci- 
entífica ou do desenvolvimento tecnológico. Pre- 
tende-se uma nota curta, realçando os resultados 
alcançados e referenciando convenientemente 
trabalhos afins, evitando a exposição longa de 
antecedentes, que deverá ser remetida para refe- 
rências. Deste modo, a originalidade que se pre- 
tende e se traduzirá pelos resultados novos que o 
autor apresenta é a marca fundamental de uma 
nota científica. 


3 Sempre que tal se ofereça conveniente, a TÉC- 
NICA ouvirá especialistas sobre a essência e es- 
trutura dos trabalhos submetidos a publicação, 
com natural excepção dos Artigos por Convite. 
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Quando o autor ou autores forem alunos é ne- 
cessário que do trabalho conste a indicação de 
um supervisor de reconhecida idoneidade cienti- 
fica e técnica que tenha comparticipado na análi- 
se critica e prévia discussão do trabalho em 
causa. 


Os Artigos por Convite e os Artigos 
Técnico-Científicos, deverão ser escritos em lin- 
gua portuguesa. As Notas Científicas poderão 
ser escritas em inglês ou francês. 


Os autores deverão atender na elaboração dos 
ss as seguintes condições: 
os Artigos por Convite não poderão exceder 
16 páginas dactilografadas, incluindo o espaço 
reservado a ilustrações; 
os Artigos Técnico-Científicos não poderão ex- 
ceder 8 páginas dactilografadas, incluindo o 
espaço reservado a ilustrações; 
- as Notas Científicas não poderão exceder 8 
páginas dactilografadas, incluindo o espaço 
reservado a ilustrações; 
nome exacto com que desejam figurar; 
- ordem de entrada dos autores, havendo mais 
do que um; 
instituição(ões) a que pertençam e que dese- 
jem referir, assim como posição que nela(s) 
ocupam, indicando claramente o nome dessa 
instituição, morada e sigla. 


Os trabalhos devem ser dactilografados a dois 
espaços, em folhas modelo A4, devidamente 
numeradas. 


Os trabalhos devem ser enviados em triplicado, 
recomendando-se aos autores que conservem 
uma cópia. 


Conjuntamente com o trabalho, devem ser envi- 
ados: 


título do trabalho em português e inglês; 
resumo em português e abstract em inglês, 
pelo menos; 

- até um máximo de seis palavras-chave, em 
português e inglês, pelo menos, que identifi- 
quem claramente a natureza do tema aborda- 
do no artigo; 

- nome, direcção e telefone do autor e demais 
dados que facilitem o processo de revisão de 
provas. 


Os autores deverão assegurar-se, antes de en- 
viarem os originais, que estes não contenham 
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erros, falta de dados, etc., pois que, na fase de 
revisão de provas, já não se admite inclusão de 
novo original e apenas se corrigirão erros de 
composição tipográfica. 


As referências bibliográficas serão indicadas no 
texto por ordem numérica, usando algarismos 
árabes entre colchetes [|, e a lista das mesmas, 
no final do artigo, obedecerá às seguintes nor- 
mas: 


Artigos de revistas: 


Nome(s) do(s) autor(es) e iniciais (ano de publi- 
cação entre parêntesis). Nome do artigo, Nome 
da revista, número do volume, página(s). 


Livros: 


1 


12 


Nome(s) do(s) autor(es) e iniciais (ano de publi- 
cação entre parêntesis). Título, editora, página(s). 


Convém distinguir claramente entre referências 
e outra bibliografia que se recomende. 

Trabalhos em que a matéria apresentada não 
esteja referenciada não são, em principio, publi- 
cados. Este é um ponto fundamental em que se 
insiste. 


Os quadros ou tabelas, numerados em romano, 
devem figurar no texto dactilografado, podendo 
ser encimados por um título conciso. 


Os gráficos devem, de preferência, ser dese- 
nhados a tinta da china, em folhas separadas e 
de modo a estarem prontos para serem repro- 
duzidos com a maior nitidez e qualidade gráfica 
possível. Quaisquer letras ou algarismos que 
neles figurem deverão ser suficientemente gran- 
des para permitir a sua redução que, em princi- 
pio, fica ao critério do editor. 


As folhas com os gráficos, figuras ou fotografias, 
deverão ser numeradas, indicando-se o trabalho 
a que dizem respeito. Além disso, no texto, indi- 
car-se-à o local para a respectiva inclusão. 


NOTA: A TÉCNICA salienta a necessidade de os au- 
tores, quando for o caso, obterem em tempo as de- 
vidas autorizações de reprodução de excertos de 
outras obras, sendo esta matéria da inteira respon- 
sabilidade do autor. 


A TÉCNICA oferecerá cinco exemplares da Revista, 
por cada trabalho publicado, independentemente do 
número de autores. 
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ASSOCIAÇÃO PARA O DESENVOLVIMENTO 
DO INSTITUTO SUPERIOR TECNICO 


ACÇÕES DE FORMAÇÃO REALIZADAS EM 1992/93 


“TÍTU LO PROGRAMA 


Formação avançada específica em técnicas de previsão 
Introdução às técnicas de Desenho Mecânico incluindo Técnicas da dispersão de poluentes no meio ambiental. Pretende-se 


TÍTULO | GRAMA 


PRO 


Desenho 


Projecto Introdução aos métodos e técnicas de análise e projecto 
em Eng. integrado por computador em Engenharia Mecânica, 
Cálculo Automático de Estruturas e Sistemas Mecânicos 
j e sua integração com sistema CAD/CAM. Utilizando : 
Integrado sofiware integrado comercial em diversos ambientes Industrial 
por informáticos. 
Computador Aplicações em projectos industriais. 


Formar técnicos de electrólise nas suas aplicações 
industriais, que para além de adquirirem uma formação 
Electrólise básica apropriada, fiquem aptos a desempenhar 
actividades profissionais de anodização, niquelagem, 
cromagem, extracção e refinação de metais e outras que 
incluem processos electroliticos, 


Mecânica 


DR mn Formação básica de electroquimica e ciência de 
fnforrndilica introdução à noventa e centíulo tati, estado materiais de modo a que os formandos fiquem aptos a 

; introdução às linguagens de programação, estrutura de ne sintetizar e caracterizar novos materiais (cerâmica, 
e Controlo informação, sistemas operativos utilitário, rede de de Materiais vidros, polímeros, compósitos) por via electroquímica, e 
Industrial comunicação, programação aplicada, automação e a montar e testar dispositivos avançados na área da 


controlo. energia óptica, análise quimica, etc. 


Electroquimica 


Formação de técnicos especializados em novas 
tecnologias experimentais. 

Ênfase especial será dado às técnicas relacionadas com o 
ultra-vácuo, a óptica aplicada, a espectroscopia, o 
controlo e a teoria da informação. 


Automação Introdução aos principios básicos da automação indus- 
trial, incluindo a teoria e a prática de análise e projecto 
de sistemas automálicos, informática industrial, robótica, 
comunicações, CIM É produção automatizada. 


Técnicas 
Exprimentais 
Avançadas 


é Controlo 
Industrial 


Descrição dos fundamentos do projecto e operação de 
pad equipamento têrmico industrial, abrangendo técnicas 
Utilização experimentais e modelação matemática. Especial atenção 
de Gás é dedicada ao controlo e à optimização das condições de | 
Natural funcionamento, tendo em vista o problema da | 
conservação da energia. 


Formação de técnicos especializados para o desempenho 
de funções em laboratórios industriais, especializados nas 
novas tecnologias empregues na automatização dos 
processos produtivos e no controlo de qualidade. 


Tecnologias de E 
| Técnicos 


de Laboratório 


Formação de técnicos indústriais (engenheiros de 
Depuração | processos, quimicos mecânicos é de ambiente) na 
e Reutilização | resolução de problemas de ambiente ligados à produção 
de Efluentes | industrial. Necessidade de actualização em novas 
Liquidos | técnicas de separação a serem utilizadas na depuração de 
águas residuais industriais e na sua incorporação na 
sinlese de processos de fabrico não poluentes. 


Formação especifica em ferramentas de escritório 
elecirônico. Pretende-se que os formandos adquiram 
conhecimentos teóricos e práticos, que viabilizem o uso 
de ferramentas de escritório electrónico de forma 
eficiente e eficaz no apoio à prática das engenharias. 


Ferramentas 
de Escritório 
Electrônico 

Industriais 


Formação especifica em programação de aplicações no 


) a Circuitos | pis dci ra FR | a seio 
Programação ambiente UNIX pretendendo-se que os formandos fiquem Integrados Formação integrada das tente as de projecto de a 
de Ambiente aptos a desenvolver aplicações em diferentes linguagens 5 integrados analógico-digitais para processamento de 


Analógicos- conversão de sinal susceptíveis de fabricação nas 


UNIX SR 
-Digitais tecnologias CMOS, BiCMOS e GaAs. 


e utilizando rotina gráfica para apoio à prática das 
engenharias. 


Formar jovens no domínio das novas técnologias da 
informação, tornando-os aptos a enfrentar uma carreira 
profissional como técnicos de informática, com 
capacidade para assumir responsabilidades de concepção 
de aplicações e sistemas. 


Técnicos 
Avançados 
de Informática 
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Mecânico normalização treino de interpretação de desenhos. Prática de Previsão que a par com uma formação teórica adequada, para 
Assistido de desenho dando-se especial ênfase no Desenho da Dispersão compreensão do fenómeno, is formando adquira prática 
por Assistido por Computador. Utilização de Programas de Poluentes de estudos de dispersão com o auxílio de métodos e 
: Comerciais em aplicações e projectos industriais. a modelos de previsão, adquirida em exercicios de 
Computador no Ambiente interesses representativos de problemas industriais. 


Jy INSTITUTO SUPERIOR TÉCNICO 


NOVAS IDEIAS PARA A 
UNIVERSIDADE 


A modernização da Universidade é uma condição essencial para que Portugal enfrente os 
desafios das próximas décadas, tanto no que respeita ao desenvolvimento económico, 
social e cultural como à sua capacidade de afirmação internacional. O novo contexto 
europeu, juntamente com a democratização do ensino e o novo relacionamento da 
Universidade com os agentes de desenvolvimento económico e social, têm alargado a 
âmbito deste problema para além da renovação de experiências pedagógicas. 


É necessário implementar novas modalidades de gestão académica, acompanhadas pela 
qualificação do nível de docência através da promoção de trabalhos de investigação de 
prestígio, e permitir novas formas de diálogo entre a Universidade e a sociedade. 


Este Ciclo de Seminários que o Conselho Directivo do IST promove em 1993 irá por certo, 
pôr em evidência algumas ideias pontuais e inovadoras que, espera-se, contribuirão para 
o esclarecimento de reforma da Universidade. 


CICLO DE SEMINÁRIOS -— 1993 


Já Realizados 
* MODERNIZAÇÃO ADMINISTRATIVA DA UNIVERSIDADE 
* ORGANIZAÇÃO E MÉTODOS EM GESTÃO DE CIÊNCIA E TECNOLOGIA 


* UNIVERSIDADE E SOCIEDADE 


A Healizar 

- INVESTIGAÇÃO E COMPETITIVIDADE INDUSTRIAL ..........ceneeeeeeeeeeeemreranaees . 24 Setembro 
* AVALIAÇÃO DAS UNIVERSIDADES ............ Ra LU EPE RERNE GER PA «. 22 Outubro 
+ CONTRATOS COMUNITÁRIOS BIRD essas ssa e aaa ca 12 Novembro 


INFORMAÇÕES E FORMULÁRIOS 


10-12H/14-16H 
Gabinete de Informação e de Apoio às Relações com o Exterior (GIRE) - IST 
Av. Rovisco Pais — Tel. 841 75 54 / 841 73 62 — Ext. 1850 — Fax. 849 92 42 — 1096 Lisboa Codex 


“revista de Engenharia 


Aqui poderia estar um anuncio seu 


Para mais informações contacte-nos pelos telefones 351-1-3559294/5 ou fax 351-1-3536682 | 
ou escreva-nos para: DOSSIER, comunicação e imagem, Lda., r. dos Acores nº84 - 1º Esq. P-1000 Lisboa 


revista de Engenharia 


Torne-se assinante da Técnica. 


4 números - 1400500 


| Para mais informações contacte-nos pelo telefone 351-1-8417410 ou fax 351-1-8499242 
| ou escreva-nos para: Técnica, revista de Engenharia, Av. Rovisco Pais nº1, P-1096 Lisboa CODEX 
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Agora, pode comprar! 

O BNU acaba de lançar uma nova forma de comprar casa. 
Com a "Habitação Segura", você só tem que constituir 
um Seguro, que lhe servirá, no final, para pagar 

a amortização do capital em dívida. 

Assim, só paga os juros do empréstimo. E aindáã poupa 
a dobrar... nos impostos. 

No BNU você encontra também todas as outras 
modalidades de Crédito à Habitação. 

Vá já ao Balcão BNU mais próximo e aconselhe-se 
sobre o crédito mais vantajoso para si. 

"HABITAÇÃO SEGURA" BNU 

Você quer? 

Você tem! 


EA GERAL DE DEPT 


BNU Banco Nacional Ultramarino 


A Telecom Portugal é a maior operadora portuguesa de telecomunicações, 
oferecendo a clientes nacionais e internacionais uma gama completa | 
de serviços. — 
Quaisquer que sejam as suas necessidades em telecomunicações, no WS. 
ambito de redes privativas ou públicas, de circuitos alugados ou de Ea 
serviços avançados, encontraremos e desenvolveremos as 
soluções mais eficientes e vantajosas. 


TELECOM PORTUGAL 
À ligação directa ao Séc. XXI 


